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Some conjectures on addition and multiplication of complex (real) numbers 

Apoloniusz Tyszka 

Abstract. We discuss conjectures related to the following two conjectures: 

(I) (see [7]) for each complex numbers xi,...,Xn there exist rationals yi,...,y„ G 
[_2n-i 2«-i] such that 

\/te{l,...,n}{x, = l^yi = l) (1) 

Wi,j,k e {I,..., n} {xi + Xj = Xk^yi + Vj = yt) (2) 

(II) (see [6], [7]) for each complex (real) numbers xi, . . . ,x„ there exist complex (real) 
numbers yi, . . . ,yn such that 

VzG {!,..., n} \yi\<2^ 

\/ie {l,...,n} {x, = l^yi = 1) 

Vz, j, A: G {1, . . . , n} {xi + Xj = Xk ^ yi + yj = yt) 

Vi, j, A; G {1, . . . , n} {xi ■ Xj = Xk ^ yi ■ yj = yt) 



(3) 

(4) 
(5) 
(6) 



Mathematics Subject Classification: 03B30, 12D99, 14P05, 15A06, 15A09 

Keywords: system of linear equations, system of polynomial equations, solution 
with minimal Euclidean norm, least-squares solution with minimal Euclidean norm 

For a positive integer n we define the set of equations Wn by 

Wn = {xi = 1: l<z<n}U {xi + Xj = x^ '■ ^ "^ i < j < n, 1 < k < n} 

Let S C Wn be a system consistent over C. Then S has a solution which consists of 
rationals belonging to [— (y^)"^^, (-\/5)"^^], see [7, Theorem 9]. Conjecture (I) states 
that S has a solution which consists of rationals belonging to [—2"^^, 2"^^]. 

Concerning Conjecture (I), estimation by 2*^"^ is the best estimation. Indeed, the 
system 

Xi = 1 

Xi + Xi = X2 
X2 + X2 = Xs 
X3 + Xs = Xa 



lij ry 1 I tt" 



n— 1 



Xi- 



has precisely one complex solution: (1, 2, 4, 8, . . . , 2" , 2" ). Concerning Conjecture (II) 

on— 2 

for n = 1 estimation by 2 can be replaced by estimation by 1. For n > 1 estimation 

on— 2 

by 2^ is the best estimation, see [7]. 

For each consistent system S C Wn there exists J C {1, . . . , n} such that the system 
SU{xi + Xi = Xi : i & J} has a unique solution (xi, . . . ,Xn), see the proof of Theorem 9 
in [7j. If any S" C Wn has a unique solution (xi, . . . , x„), then by Cramer's rule each Xi 
is a quotient of two determinants. Since these determinants have entries among —1, 0, 
1, 2, each Xi is rational. 

For proving Conjecture (I), without loss of generality we can assume that the 
equation xi = 1 belongs to S and all equations Xi = 1 {i > 1) do not belong to S. 
Indeed, ii i > 1 and the equation Xj = 1 belongs to S, then we replace Xj by Xi in all 
equations belonging to S. In this way the problem reduces to the same problem with 
a smaller number of variables, for details see the text after Conjecture 3. Therefore, 
for proving Conjecture (I) it is sufficient to consider only these systems S C Wn of n 
equations which satisfy the following conditions: 

S contains the equation Xi = 1 and n — 1 equations of the form x, + Xj = x^ 
{i,j,k e {l,...,n}), 

joining the equation Xi = and the aforementioned n — 1 equations of the form 
Xi + Xj = Xfc, we get a system of n linearly independent equations. 

By the observations from the last two paragraphs, the following code in MuPAD 
yields a probabilistic confirmation of Conjecture (I). The value of n is set, for example, 
to 5. The number of iterations is set, for example, to 1000. We use the algorithm which 
terminates with probability 1. For another algorithm, implemented in Mathematica, 
see [4]. 

SEED:=time() : 

r :=random(l . .5) : 

idmatrix : =matrix : : identity (5) : 

u:=linalg: :row(idmatrix,i) $i=l. .5: 

max_norm:=l : 

for k from 1 to 1000 do 

a:=linalg: :row(idmatrix, 1) : 

rank:=l : 

while rank<5 do 

m:=matrix(u[r()] )+matrix(u[r()] ) -matrix (u[r()] ) : 



al:=linalg: : stackMatrix(a,m) : 

rankl :=linalg: :rank(al) : 

if rankl > rank then a:=linalg: :stackMatrix(a,m) end_if: 

rank:=linalg: : rank (a) : 

end_while : 

x: = (a"-l)*linalg: : coKidmatrix, 1) : 

max_norm : =max (max_noriii , norm (x) ) : 

print (max_noriii) : 

end_f or : 

Conjecture (I) holds true for each n < 4. It follows from the following Observation 1. 

Observation 1 ([7, p. 23]). If n < 4 and {xi, . . . , x„) G C" solves 5", then some (fl, . . . , x^ 
solves 5", where each Xi is suitably chosen from {xi, 0, 1, 2, |} fl {r G Q : \r\ < 2"~^}. 

Let Ax = b be the matrix representation of the system S, and let A^" denote 
Moore-Penrose pseudoinverse of A. The system S has a unique solution xq with 
minimal Euclidean norm, and this element is given by xq = A'I'b, see [5, p. 423]. 

For any system S C Wn, a vector x G C" is said to be a least-squares solution if x 
minimizes the Euclidean norm of Ax — b, see [1, p. 104]. It is known that xq = A^b 
is a unique least-squares solution with minimal Euclidean norm, see [1, p. 109]. 

Since A has rational entries (the entries are among —1, 0, 1, 2), A''" has also rational 
entries, see [2, p. 69] and [3, p. 193]. Since b has rational entries (the entries are among 
and 1), xo = A^b consists of rationals. 

Conjecture 1. The solution (The least-squares solution) xq consists of numbers be- 
longing to [-2"~\2"-i]. 

Conjecture 1 restricted to the case when card S < n implies Conjecture (I). The 
following code in MuPAD yields a probabilistic confirmation of Conjecture 1 restricted 
to the case when card S < n. The value of n is set, for example, to 5. The number of 
iterations is set, for example, to 1000. 

SEED:=time() : 

r :=randoiii(l . .5) : 

idmatrix : =matrix : : identity (5) : 

u:=linalg: :row(idmatrix,i) $i=l. .5: 

max_norni:=l : 

for k from 1 to 1000 do 



b: = [l]: 

c:=linalg: :row(idmatrix, 1) : 

for w from 1 to 4 do 

h:=0: 

hl:=r() 

h2:=r() 

h3:=r() 

ml : =matrix (u [hi] ) 

m2 : =matrix (u [h2] ) 

m3 : =matrix (u [h3] ) 

m:=ml+m2-m3: 

c:=linalg: :stackMatrix(c,m) : 

if h3=h2 then h:=l end_if : 

b:=append(b,h) : 

a:=linalg: :pseudoInverse(c) : 

x:=a*matrix(b) : 

max_norm : =max (max_norm , norm (x) ) : 

end_f or : 

print (max_norm) : 

end_f or : 

The following Conjecture 2 implies Conjecture (I), see [7J. 

Conjecture 2 ([3)- Let B be a matrix with n — 1 rows and n columns, n> 2. Assume 
that each row of B, after deleting all zeros, forms a sequence belonging to 

{(1), (-1,2), (2,-1), (-1,1,1), (1,-1,1), (1,1,-1)} 
We conjecture that after deleting any column of B we get the matrix whose determinant 
has absolute value less than or equal to 2"~^. 

Conjecture 3. If a system S C Wn has a unique solution (xi, . . . , Xn), then this solu- 
tion consists of rationals whose nominators and denominators belong to [—2"^^, 2"^^]. 

Conjecture 3 implies Conjecture (I). The MuPAD code below confirms Conjecture 3 
probabilistically. As previously, the value of n is set to 5, the number of iterations is 
set to 1000. We declare that 

{z G {1, 2, 3, 4, 5} : the equation Xj = 1 belongs to 5} = {1}, 

but this does not decrease the generality. Indeed, (0, . . . , 0) solves 5* if all equations 
Xi = 1 do not belong to S. In other cases, let 



I = {i E {I, . . . ,n} : the equation Xj = 1 belongs to 5"}, 

and let i = min(/). For each j E I we replace Xj by Xi in all equations belonging to S. 
We obtain an equivalent system S with n — card(/) + 1 variables. The system S has a 
unique solution (fi, . . . , tn-card{/)+i), and the equation Xj = 1 belongs to S if and only 
if j = i. By permuting variables, we may assume that i = 1. 

SEED:=time() : 

r :=random(l . .5) : 

idmatrix : =matrix : : identity (5) : 

u:=linalg: :row(idmatrix,i) $i=l. .5: 

abs_numer_denoiii: = [l] : 

for k from 1 to 1000 do 

c:=linalg: :row(idmatrix, 1) : 

rank:=l : 

while rank<5 do 

m:=matrix(u[r()] )+matrix(u[r()] ) -matrix (u[r()] ) : 

if linalg: :rank(linalg: : stackMatrix(c,m))>rank 

then c :=linalg: :stackMatrix(c,m) end_if: 

rank:=linalg: :rank(c) : 

end_while : 

a:=(c"-l)*linalg: : col (idmatrix, 1) : 

for n from 2 to 5 do 

abs_numer_denom:=append(abs_numer_denom,abs(numer(a[n] ) )) : 

abs_numer_denom:=append(abs_numer_denom,abs(denom(a[n] ) )) : 

end_f or : 

abs_numer_denom:=listlib: :removeDuplicates(abs_numer_denom) : 

print (max (abs_numer_denom)) : 

end_f or : 

The MuPAD code below completely confirms Conjecture 3 for n = 5. We declare 
that 

{i E {1, 2, 3, 4, 5} : the equation Xj = 1 belongs to S} = {1}, 
but this does not decrease the generality. 

p: = []: 

idmatrix : =matrix : : identity (5) : 

u:=linalg: :row(idmatrix,i) $i=l. .5: 



for rl from 1 to 5 do 

for r2 from 1 to 5 do 

for r3 from 1 to 5 do 

ml : =matrix (u [rl] ) : 

m2 : =matrix (u [r2] ) : 

m3 : =matrix (u [r3] ) : 

m:=ml+m2-m3: 

p:=append(p,m) : 

end_f or : 

end_f or : 

end_f or : 

p:=listlib: :removeDuplicates(p) : 

p:=listlib: : setDifference(p, [linalg: :row(idmatrix, 1)] ) 

abs_numer_denom:=[] : 

si :=nops(p)-l : 

s2:=nops(p)-2: 

s3:=nops(p)-3: 

for n3 from 1 to s3 do 

w2:=n3+l: 

for n2 from w2 to s2 do 

wl:=n2+l: 

for nl from wl to si do 

wO:=nl+l: 

for nO from wO to nops(p) do 

c3:=linalg: : stackMatrix (linalg: :row(idmatrix, 1) ,p[n3]) : 

c2:=linalg: : stackMatrix (c3,p[n2] ) : 

cl:=linalg: : stackMatrix (c2,p[nl] ) : 

c :=linalg: :stackMatrix(cl,p[nO] ) : 

if linalg: :rank(c)=5 then 

a: = (c"-l)*linalg: : coKidmatrix, 1) : 

for n from 2 to 5 do 

abs_numer_denom:=append(abs_numer_denom,abs(numer(a[n] ) )) : 

abs_numer_denom:=append(abs_numer_denom,abs(denom(a[n] ) )) : 

end_f or : 

abs_numer_denom : =listlib : : removeDuplicates (abs_numer_denom) 



end_if : 

end_f or : 

end_f or : 

abs_numer_denoin : =sort (abs_numer_denoiii) : 

print (abs_numer_denom) : 

end_f or : 

end_f or : 

The following Conjecture 4 implies Conjecture (I), because each consistent system 
S C Wn can be enlarged to a system S C Wn with a unique solution (xi, . . . , Xn) and 
(a:i,...,x„)eQ". 

Conjecture 4. Let rationals xi, . . . ,x„ satisfy |xi| < \x2\ < • • • < |a;„|, and for each 

yi,...,yneQ, if 

\/ie{l,...,n} {xi = l^yi = 1) 
and 

Vi, j, A; G {1, . . . , n} {xi + Xj = Xk ^ Vi + Vj = yk), 
then (xi, . . . , Xn) = {yi, • • • , yn)- We conjecture that for each i G {1, . . . , n — 1} the 
inequality \xi\ > 1 implies Ixj+il < 2 ■ |xj|. 

Concerning Conjecture 4, the assumption |xj| > 1 is necessary to state that |xj+i| < 
2 ■ |xj|. As a trivial counterexample we have (0, 1), (|, |, 1, |, 2, 3) is a counterexample 
which consists of positive rationals alone. The MuPAD code below probabilistically 
confirms Conjecture 4. The value of n is set, for example, to 5. The number of iterations 
is set, for example, to 1000. 

SEED:=time() : 

r :=randoiii(l . .5) : 

idmatrix : =matrix : : identity (5) : 

u:=linalg: :row(idmatrix,i) $i=l. .5: 

max_ratio :=1 : 

for k from 1 to 1000 do 

a:=linalg: :row(idmatrix, 1) : 

rank:=l : 

while rank<5 do 

m:=matrix(u[r()] )+matrix(u[r()] ) -matrix (u[r()] ) : 

al:=linalg: : stackMatrix(a,m) : 

rankl :=linalg: :rank(al) : 



if rankl > rank then a:=linalg: :stackMatrix(a,m) end_if: 

rank:=linalg: : rank (a) : 

end_while : 

x: = (a"-l)*linalg: : coKidmatrix, 1) : 

xx:=[max(l,abs(x[i])) $i=1..5]: 

xxx:=sort(xx) : 

maxratio:=max( [xxx[i+l] /xxx[i] $i=l . .4] ) : 

max_ratio :=max(max_ratio,maxratio) : 

print (max_ratio) : 

end_f or : 

For a positive integer n we define the set of equations En by 

En = {xi = l: 1 <i<n} U 
{xi + Xj = Xk '■ ^ ^ i < j < n, 1 < A; < n} U {xi ■ Xj = x^ '■ 1 < ^ < i < "^j 1 <k <n] 

Let T C i?„ be a system consistent over C (over R). Conjecture (II) states that T 
has a complex (real) solution which consists of numbers whose absolute values belong 

on— 2 

to [0,2 ]. Both for complex and real case, we conjecture that each solution of T 
with minimal Euclidean norm consists of numbers whose absolute values belong to 

on — 2 

[0, 2 ]. This conjecture implies Conjecture (II). Conjecture (II) holds true for each 
n < 4. It follows from the following Observation 2. 

Observation 2 ([7, p. 7]). If n < 4 and (xi,...,x„) G C" (M") solves T, then 
some (fi,...,5^) solves T, where each Xi is suitably chosen from {xj, 0, 1, 2, |}n 

on— 2 

{zeC (M) : 1^1 <2^ }. 

Let us consider the following four conjectures, analogical conjectures seem to be 
true for R. 

(5a) If a system 5" C i?„ is consistent over C and maximal with respect to inclusion, 
then each solution of 5" belongs to 

, , on— 2 , , on— 2, 

{(a;i,...,x„) G C" : |xi| < 2^ A ... A |a;„| < 2^ }. 



(5b) If a system 5" C £"„ is consistent over C and maximal with respect to inclusion, 
then S has a finite number of solutions (xi, . . . , a;„). 

(5c) If the equation Xi = 1 belongs to 5" C i?„ and S has a finite number of complex 
solutions (xi, . . . , x„), then each such solution belongs to 

, , on— 2 , , on— 2, 

{(a;i,...,x„) G C" : |xi| < 2^ A ... A |x„| < 2^ }. 



(5d) If a system S <^ En has a finite number of complex solutions (xi, . . . , Xn), then 
each such solution belongs to 



{(xi,...,x„)eC": \x^\ <2^ 



in— 1 



A ... A |x„l < 2" 



in— 1 , 



}• 



Conjecture 5a strengthens Conjecture (II) for C. The conjunction of Conjectures 
5b and 5c implies Conjecture 5a. 



Concerning Conjecture 5d, for n = 1 estimation by 2 

in— 1 



in— 1 



estimation by 1. 


For n > 1 estimation by 


2^ 


the system 


• 








Xi + Xi 


= X2 






Xi ■ Xi 


= X2 




< 


X2 ■X2 


= X3 






X3-X3 


= X4 






Xn—1 ' Xn—1 


^ Xn 



can be replaced by 
is the best estimation. Indeed, 



on— 2 on— 1 

has precisely two complex solutions, (0, . . . , 0), (2,4, 16, 256, . . . , 2 ,2 ). 

The following code in MuPAD yields a probabilistic confirmation of Conjectures 5b 
and 5c. The value of n is set, for example, to 5. The number of iterations is set, for 
example, to 1000. 

SEED:=time() : 

p:=[v-l,x-l,y-l,z-l] : 

var :=[l,v,x,y,z] : 

for i from 1 to 5 do 

for j from i to 5 do 

for k from 1 to 5 do 

p : =append (p , var [i] +var [ j ] -var [k] ) : 

p : =append (p , var [i] *var [ j ] -var [k] ) : 

end_f or : 

end_f or : 

end_f or : 

p:=listlib: :removeDuplicates(p) : 

max_abs_value : =1 : 

for r from 1 to 1000 do 

q:=combinat: : permutations : : random (p) : 

syst : = [t-v-x-y-z] : 
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w:=l : 

repeat 

if groebner: : dimension (append (syst ,q[w] ))>-! 

then syst :=append (syst ,q[w] ) end_if: 

w:=w+l : 

until (groebner : :dimension(syst)=0 or w>nops(q)) end: 

d:=groebner: : dimension (syst) : 

if d>0 then print ("Conjecture 5b is false") end_if : 

if d=0 then 

sol : =numeric : : solve (syst) : 

for m from 1 to nops(sol) do 

for n from 2 to 5 do 

max_abs_value :=max(max_abs_value,abs(sol[m] [n] [2])) : 

end_f or : 

end_f or : 

end_if : 

print (max_abs_value) ; 

end_f or : 

If we replace 

p: = [v-l,x-l,y-l,z-l] : by p: = []: 

var : = [1 , V , X , y , z] : by var : = [u , v , x , y , z] : 

max_abs_value :=1 : by max_abs_value:=0: 

syst :=[t-v-x-y-z] : by syst :=[t-u-v-x-y-z] : 

for n from 2 to 5 do by for n from 2 to 6 do 

then we get a code for a probabilistic confirmation of Conjecture 5d. 

It seems that for each integers Xi, . . . ,Xn there exist integers yi, . . . , y„ G [— 2"~\ 2"~^] 
with properties (1) and (2), cf. [7, Theorem 10]. However, not for each integers 
xi, . . . ,Xn there exist integers yi, . . . ,yn with properties (3)-(6), see [7, pp. 15-16]. 

The author used MuPAD Pro 4.0.6. SciFace Software GmbH & Co. KG, the 
maker of MuPAD Pro, has been acquired by The MathWorks, the maker of MATLAB 
technical computing software. MuPAD Pro is no longer sold as a standalone product. 
A new package, Symbolic Math Toolbox 5.1 requires MATLAB and provides large 
compatibility with existing MuPAD Pro applications. 
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